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Abstract 

Both small world models of random networks with occasional long range connec- 
tions and gossip processes with occasional long range transmission of information 
have similar characteristic behaviour. The long range elements appreciably reduce 
the effective distances, measured in space or in time, between pairs of typical points. 
In this paper, we show that their common behaviour can be interpreted as a product 
of the locally branching nature of the models. In particular, it is shown that both 
typical distances between points and the proportion of space that can be reached 
within a given distance or time can be approximated by formulae involving the limit 
random variable of the branching process. 
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1 Introduction 

Moore & Newman [7] introduced a continuous analogue of the Watts & Strogatz [8] 
"small world" model. In this model, a random number of chords, with Poisson distribu- 
tion Po(Lp/2), are uniformly and independently superimposed as shortcuts on a circle 
of circumference L. Distance is measured as usual along the circumference, but chords 
are deemed to be of length zero, and interest centres on finding the statistics of shortest 
path distances between pairs of points. A closely related model, the "great circle model", 
was introduced somewhat earlier by Ball et al. [2], in the context of epidemics; here, dis- 
tance between points translates into time taken for one infected person to infect another. 
In Barbour & Reinert pE], assuming the expected number Lp/2 of shortcuts to be large, we 
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proved a distributional approximation for the distance between a randomly chosen pair 
of points P and P', and gave a bound on the order of the error, in terms of total variation 
distance. We also showed that analogous results could be proved in higher dimensions 
by much the same method, when the circle is replaced by a sphere or a torus. It turns 
out that the reduction in the typical distances between pairs of points that results from 
introducing shortcuts is still substantial, but less dramatic than in one dimension. 

More recently, Chatterjee & Durrett [5] studied a model for the spread of gossip that 
is the continuous analogue of one of a number of models discussed in Aldous p]. Here, 
information spreads locally from an individual to his neighbours on the two-dimensional 
torus, and also occasionally to other, randomly chosen members of the community. Thus 
a disc of informed individuals, centred on an initial informant, grows steadily in the torus; 
long range transmissions of information occur in a Poisson process, whose rate is propor- 
tional to the area (number) of informed individuals, and any such transmission contacts 
a randomly chosen point of the torus, initiating a new disc of informed individuals. The 
distinction between this model and the corresponding two dimensional model in j3] is 
that, in the gossip model, the Poisson process runs at a rate proportional to the area of 
the currently informed region; in [3], where the Poisson number of shortcuts is consid- 
ered to be fixed in advance, the Poisson process corresponds to a process of discovery of 
shortcuts, and its rate is thus proportional to the length of the boundary of the informed 
region. 

Here, we consider the development of such a process on a smooth closed homogeneous 
Riemannian manifold C of dimension d, such as a sphere or a torus, having large finite 
volume \C\ =: L with respect to its intrinsic metric. We then assume that, around each 
point P of C, there is a collection of closed subsets K,(P,s), s > 0, that are balls of radius s 
with respect to a metric dc that makes C a geodesic space, and with (intrinsic) volumes 
v s (K.) := \fC(P,s)\ ~ s d v{JC) as s — > 0, for some v(]C) > 0; s is thought of as time, and 
{w(/C)} 1//rf as a (linear) speed of propagation. We shall in particular assume that 

\s- d v a {K)/v{K)-l\ < c g (s{v(lC)/L} 1 / d y*, s>0, (1.1) 

for some 7 g > 0. The metric dc need not be the same as the intrinsic metric; for instance, 
on the torus, we could consider rectangular as well as circular neighbourhoods. Two such 
subsets K,(P,t) and JC(Q,u) then intersect when P e JC(Q,t + u), or, equivalently, when 
Q G JC(P, t+u), so that the probability of intersection if P and Q are chosen independently 
and uniformly distributed on C (with respect to the intrinsic volume), is given by 

q L {t,u) := L~ l v t+u {1C) = L-yjCXt + iO^l + ifcftu)), (1-2) 

where 

\R L {t,u)\ < c g {{t + u){v{)C)/L) 1 / d y°. (1.3) 

The set JC(P,s) denotes the set of points 'locally' contacted after time s has elapsed 
following an initial 'long range' contact at P, thought of as 'islands' in C, and the complete 
set of contacts Yp (t) at time t is the union of these sets growing from an initial point P®, 
and from all long range contacts made before t. The rate at which long range contacts are 
made is proportional either to the area of the boundary of Y(t) (small world processes) 
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or to its volume (gossip processes), and we denote the constant of proportionality by p; 
long range contacts are made to independently and uniformly chosen points of C. 

It is clear that, at least for a while, such a process can be closely approximated using a 
Markovian growth and branching process X* taking values in U,->iC J ', where C denotes the 
closed subsets of C. Long range contacts are made to independently and uniformly chosen 
points of C, and, at time s after a contact at P*, the set /C(P*, s) is one component of the 
state of the process. Thus the state of the branching process at time t can be encoded in 
terms of the position P * e C of the initial individual and the number n(t) of contacts that 
have taken place up to time t, together with the pairs (r*, P*), . . . , (t*^, P*( t )) € R+ x C 
of times and positions of the contacts, giving 

X*(t) = (J*(0,t),J*(l,t),...,r(n(t),t)), 

where J*(j,t) : = /C(P*,t - r*) if t* < t and J*(j,t) = otherwise; we take r * = 0. 
In 'gossip' models, new contacts are made at a rate proportional to the current volume, 
which, for the process X* is given by 

n(t) n(t) 

V*(t) := J>_ r ;(/C) ~ J2(t-T*) d v(>Cy, 

3=0 3=0 

in 'small world' models, the rate is proportional to the derivative of the volume. 

If Yp(t) C C denotes the state of the actual process of interest, then it can initially 
be approximated by taking P * = P, and then forming the union 

n(t) 

Y k{t) ■= \Jj*U,t). (1.4) 

3=0 

However, as soon as the union is not disjoint, the actual rate of contacts decreases, and the 
two processes diverge. To be more precise about the definition of Y — Yp* in terms of X*, 
we augment the encoding of X* by adding to each pair (r*, P*) both the position Qj e C 
and the index Ij of the island from which the j-th contact was made, together with a 
mark Gj which is if the j-th island is 'real', and 1 if it is a 'ghost' island. Gj is specified 
by setting Gj = 1 (island j is a ghost) if G tj = 1, or if P* e Y(r*—), or if Qj ^ J Ijtj , 
where 

J u ■= {P eC: d c (P, P*) < d c (P, P£) for all m < j such that G m = 0}. 

To explain these restrictions, note that new branching process contacts (designated by 
Gj = 0) contribute to Y only if they are originate from a real island, and start at a point 
outside Y, and multiple counting from overlapping islands is also to be prohibited. Thus 

Y(t) := |J P(j,t), (1.5) 

where Jy(t) '■— {j > 0: r* < t, Gj = 0}, and its volume V(t) and area A(t) are in this 
way functions of {(T*,P*,Qj,Ij): j > 0,r* < t}. In view of this construction, the two 
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processes Y and Y* are coupled so as to make the identification of Y(t) and Y*(t) exact, 
until the union in (jl.4p ceases to be disjoint. 

In [I], the distribution of inter-point distances is determined by running two such 
branching processes from randomly chosen initial points, each for a length of time t* at 
which the mean number of overlaps in (jl.4p is of order 0(1). At this time, conditionally 
on the contact times in the two branching processes, the number of permissible overlaps 
between the sets (11.41) — cases in which an island in one branching process is contained 
within an island in the other could not have arisen in the actual process — has ap- 
proximately a Poisson distribution, and the distance between the initial points is greater 
than 2t* if there are no permissible overlaps. In this way, and by varying the choice of t* 
appropriately, the distribution of inter-point distances can be approximated, without ever 
having to go into the dependence structure that becomes important in the process Y at 
times significantly larger than t*. In contrast, Chatterjee & Durrett [5J go beyond the 
branching phase in the analysis of Y in their two-dimensional gossip model, and are able 
to prove a conditional law of large numbers for the fraction of the torus contained in Y, 
given the outcome of the branching phase. They also establish the asymptotics of the 
first time at which Y is the whole torus. 

The main result of this paper is Theorem I3.2[ which establishes a conditional law of 
large numbers, together with a rough error bound, for the time evolution of the covered 
fraction of C, in the setting of a quite general small world or gossip process. In the 
particular case of the torus, this extends the limit law proved by Chatterjee & Durrett [5J, 
by providing an estimate of the approximation error that is uniform for all time. Our 
argument, as in [4j, avoids any detailed analysis of Y beyond the branching phase. We 
first note that the asymptotics of the mean fraction KV(t)/L of C that is contained 
in Y, for t in the relevant range (corresponding to 2t* above), can be deduced for small 
world processes from Theorem 4.2 of [I], and we show that an analogous theorem holds 
also for gossip processes. The same argument also gives an asymptotic expression for 
E{V(i)/L | J^s}, the conditional expectation of the covered fraction at t, given the history 
of the process up to time s. The law of large numbers is then proved by using an argument 
of much the same flavour, now involving three independent branching processes in their 
early phases, to derive an asymptotic formula for E{[V(t)/L] 2 | J 7 ,}. This is shown to be 
asymptotically the same as [E{V(t)/L | ^s}] 2 , for large enough s, from which it follows 
that the conditional variance of V(t)/L, given the information in jF s , is small, and hence 
that the value of V(t)/L is (almost) fixed. An analogous argument is used, for instance, in 
Ball, Sirl & Trapman |3J, where they show that the proportion of susceptibles infected an 
epidemic in a large population is close to a fixed value, provided that the epidemic is a large 
one. A by-product of our argument is to identify the solution h to a particular integral 
equation, that appears in Aldous pQ and also plays a substantial part in the formula given 
by Chatterjee Sz Durrett [5J, in terms of the Laplace transform of the branching process 
limit random variable W. The paper concludes by extending the results to finite subsets 
of homogeneous manifolds, such as rectangles in R 2 . 
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2 The branching phase 

As in [3], we base our analysis of the coverage process on the pure growth Markov branch- 
ing process X*, which has neighbourhoods with centres independently and uniformly 
positioned in C. In this section, to describe the behaviour of such processes, we specialize 
to the case of 'flat' manifolds, such as tori, in which 

v.{K) = s d v{K), s>0, (2.1) 

so that the constant c g in (jl.ip is zero. We later show that processes in which c g cannot 
be taken to be zero can be dealt with by bounding them between processes satisfying 
condition (12. ip that are close enough for our purposes. 

We begin by defining M (t) := max{j > 0: r* < t} to be the number of contacts in 
the branching process up to time t, and 

M (t) 
j=0 

to be the sum of the Z'th powers of their 'radii'. In the small world process, 

M (t) = M o (0) + Z (pv(K)d J^M d ^{u)d^ = M o (0) + Z{pv{K)[M d (t) - M d (0)\), 

(2.3) 

where Z is a unit rate Poisson process, and in the gossip process 

M (t) = M o {0)+Z (pv(K.) J M d {u)d^j = M (0)+Z(p(rf+l)- 1 t;(/C)[M d+1 (t)-M d+1 (0)]); 

(2.4) 

in both cases, the intensity p may depend on L. The remaining evolution is governed by 
the differential equations 

— Mi(t) = M (t) for a.e. t, 

d (2 ' 5) 
j t Mi{t) = iM^it), i > 2. 

These equations can be rewritten in clearer form by defining Hi(t) := Mi(t)X l /i\, for A 
to be suitably chosen, in which case (I2.5P reduces to 

-fHtit) = XH (t) for a.e. t, 

dt d (2.6) 
—Hi(t) = XH^it), %>1- 

for the small world process, we have 
H (t) = M Q (0) + Z(d\pv()C)\- d [H d (t)-H d (0)}) = M o (0) + Z(H d (t) - H d (0)), (2.7) 
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if A = A := (<i!pt>(/C)) 1 / d , and, for the gossip process, we have 

H (t) = M (0) + Z(rf!pt;(/C)A- d - 1 [iJ rf+1 (t)-i7 rf+1 (0)]) = M o (0) + Z(H d+1 (t) -H d+1 (0)), 

(2.8) 

if A = Ao := (rf!pt>(/C)) 1 /^ +1 \ Note that, since p may depend on L, so may Ao- 
Remark 2.1 The time-scaled process H(u) := H(u/X) actually satisfies 

^Hi{t) = H^t), z>l; H (t) = M (0) + Z(H r{d) (t)-H r{d) (0)), (2.9) 

where r(d) = d for the small world and d + 1 for the gossip process. Thus, apart from 
a time change, the processes are the same for all A. Despite this, we retain Ao in the 
subsequent discussion, in order to emphasize the connection with the original process. 

In either case, the equations for H = (Hi, H2, . . . , H r ) T are of the form 

^ = \ C r (H + he r ) = X C r [I + {h/H r )e r {e r ) T )H, (2.10) 

where r = r(d), C r is the r-dimensional cyclic permutation matrix, e % denotes the z-th 
coordinate vector, and 

h(t) := H (t)-H r (t) = Z(H r (t)-H r (0))-H r (t) + M (0). (2.11) 

Without the perturbation h, H would have asymptotically exponential growth at rate Ao, 
and the ratios of its components would all tend to unity, since the dominant eigenvalue 1 
of C r corresponds to the right eigenvector 1. For the arguments to come, it will be 
important to show that, with high enough probability, the asymptotic effect of h is just 
to multiply H by some random constant — a branching random variable W — which is 
not too big. Unless otherwise specified, we henceforth take Mo(0) = 1 and M/(0) = for 
all / > 1, so that we start with just one point Pq at t = 0. 

2.1 Growth bounds for the branching process 

Using the maximum norm || • || for r(d)-vectors, it follows immediately from (I2.10p that 

j t \\H(t)\\ < Xou(t)\\H(t)\\, 
with u(t) := {1 + (h(t)/ H r ^(t)) + }, so that, by a Gronwall argument, 

\\H(t)\\ < ||#(t )||exp jAojT u{v)d v y (2.12) 

for any < to < t. Thus, in order to bound the growth of H, we shall need to control 
the quantity h(t)/ H r ^(t), which is itself a function of the Poisson process Z. To do so, 
we begin with the following lemma, which controls the extreme fluctuations of Z. 
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Lemma 2.2 Let Z be a unit rate Poisson process. Then we have the following bounds, 
uniformly in t > 1 : 



(1) P 

(2) P 

(3) P 



supw l Z{u) > 2 

u>t 



< C\ e 



sv^u l/z \u~ l Z(u) - 1| > 4 

u>t 



-t/lA. 



< c 2 e 



inf u~ L Z(u) < 1/2 



u>t 



< c 3 e 



-t/32 



Furthermore, for any U > 1 and < i] < 1/3 suc/j that U(2 V — 1) > 421og2 ; 



(4) P 



sup(w V 1) 

u>0 



4(1-H?) 



|Z(u) -u| > U 



< c 4 e" 



-(7/28 



/or a constant C4. 

Proof: For any £, e > 0, set Uj := £(1 + e) 3 ', j > 0. Then it is immediate that 

sup u~ x Z{v) < Z(v,j+i)/uj. 

Hence, and by the Chernoff inequalities ([6], Theorem 3.2), 

P sup u ^Z(u) > 1 + 2e\ < Pfuji^u^l + e)) > 1 + 2e] 

< exp{-e 2 M j /(2 + 3£)} 

= exp{-e 2 M i _i/(2 + 3e)} exp{-£ 3 «,_i/(2 + 3e)}, 



and 



P 



inf u~ x Z(u) < 1 - 2e < P[{(1 + e)^} -1 ^^) < 1 - 2e] 

< exp{-e 2 u i /2} 

= exp{— e 2 Uj-i/2} exp{— e 3 Uj_i/2}. 



Adding over j > 0, it thus follows that 



P 



sup|u _i Z(V) - 1| > 2e 



u>t 



< C(e,t)exp{-e 2 t/(2 + 3e)}, 



(2.13) 



with C(e,t) := 2/{l — e - £3 */( 2 + 3£ )}. Taking e = 1/2 gives the first inequality, with 
Ci := C(|,l); taking £ = 1/4 gives the third, with c 3 = (7(4,1). For the second, with 
t > 1, e = i -1 / 3 gives, in particular, 



P 



sup u l ^\u' l Z{u) - 1| > 4 

t<u<St 



< (7(1, 1) exp{-t 1 /3/5} j 
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and thus 



P 



supn 1/3 |n _1 Z(n) - 1| > 4 



u>t 



< C(l,l)^exp{-2^ 1 / 3 /5} < csexpl-t 1 / 3 ^}, 



with c 2 := C(l, 1)/(1 - e" 1 / 5 ). The fourth inequality is a little trickier. Taking t > 1 and 

£ = f//{2(2t)^ 1 -^}, we have 



P 



sup u^-'rtlu-^u) - 1| > U 



t<u<2t 



< C(e,t) exp{-2- 3+r 'f/ 2 tV(2 + 3Z7/2)} 

< C(£,t)exp{-C/^/28}. 



For this choice of £, e 2 t increases with t, but e 3 t decreases; however, it is not difficult to 
show that 



C(e,t) < C'(e,t) := lOe Q V -L) 



5 e 3 t ; ' 

uniformly in t, U > 1. Set g(f) := C"(e, t) exp{-C/^/28}. Then, in the sum Ej>o<?( 2i )> 
the ratios of successive terms are at most 

2 (1 - 3r?)/2 exp{-f/(2 r ' - l)/28} < v / 2exp{-f/(2 T ' - l)/28} < 1/2, 

by assumption, so that Ej>cW( 2j ) ^ 2( l( l ) ^ 2 11 / 2 lOe exp{-[//28}. Since P[Z(1) > 
U — 1] < ce -C/ / 14 by an exponential moment inequality, with c := xe 11-1 and x = e 1 / 14 , 
the proof of the fourth inequality is complete. ■ 

Based on this lemma, we can now prove growth bounds for the Markov branching 
process. Here, we allow for quite general initial conditions. For ease of reference, for any 
K > 1 and < rj < 1, we define the events 

A ( t ■= {e- x n\H(s)\\<K}; A%. := {H r(d) (s) > K}; (2.14) 
A%\ s := \ sup (mV1)-^|Z( W )- M |<^ 1 -"» ; (2.15) 

[0<u<H r(d) (s) J 

A'(K,s) := {esxp{-Xo{t-s)(l + e K )}\\H(t)\\ < \\H(s)\\ for alU > s}, (2.16) 
where Ek := 5i^ _1//3 , and we write 

Ak, s ■■= 4l s nA% s nA% eK;S e T„ (2.17) 

where T s denotes the history of X* up to time s, and 9 := Cae 1 ^ 80 , with C a as defined 
below. 

Theorem 2.1 For any K > 1 and any < s < t, we nave 

(1) P[exp{-Ao(t-s)(l + £A-)}||^(*)|| < max{C a K, ||#(s)||} for all t > s] 

for suitable constants C a '■= 3 exp{(r(rf)!) 1 / r ^} and c a . Furthermore, 

(2) F[A'(K,s)\ F 8 nA%\] > l-c 2 e- KirA '\ 
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Proof: For K > 0, define tk, s '■= inf{t > s: H r /^(t) > K}, and suppose first that 
tk, s > s. Then H r ^){t) < K if s < t < tk, s , and thus, for all such t, 

H (t) < 1 + Z(K) < 3K, 

by Lemma I2T27 1 ) . on a set A\{K) of probability at least 1 — c\ exp{— K/1A}. Hence, since 
from the definition of Hi(t) and by Holder's inequality, we have 

Hi(t) < H r[d) (tf r ^ H (ty-^ d \ l<i<r(d), (2.18) 

it follows that, on A\{K\ 

\\H(t)\\ < 3exp{(r(d)!) 1 / r W} J K' for all s < t < r K>s . (2.19) 

Now, from Lemma [2.2( 2). it follows that 

u(t) = l + (h(t)/H r{d) (t))+ < max{l,K- x + Z(H r(d) (t))/H rW (t)} 

< 1 + 5K- 1/3 = l + e K 

for all t > tk :S7 on an event A 2 (K) of probability at least 1 — C2e~ Kl/3 ^ 5 . By (I2.12p . this 
implies that, on A 2 (K) fl Ai(K), for all t > t KjS , 

\\H(t)\\ < \\H(T KtS )\\ex V {\ (t-s)(l + e K )}. (2.20) 

If tk, s > s, by (12.191) . this in turn implies that, on A 2 (K) fl Ai(K), 

\\H(t)\\ < 3exp{(r(rf)!) 1 /^}irexp{A (t- S )(l + £x)} (2.21) 

for all t > s; if t KjS = s, we simply have < ||if(s)|| exp{A (t — s)(l + Ek)}- This 

establishes Part 1. 

For Part 2, if H r u)( s ) > K, it follows as above that 

F[u(t) < l + e K ior allt> s\ T s ] > 1 - c 2 e~ xl/3/5 , 
and, if this is the case, then 

exp{-A (t-s)(l + £ K )}||F(£)|| < \\H(s)\\ for all t > s 
follows from fl2TT2j) . ■ 

Corollary 2.3 Given any e > 0, there exists a random variable H e such that 

(1) \h(t)\ < # e exp{|A (l + e)t} a.s. for all t > 0. 
In addition, for any K > 1, 

(2) E{fh(t)\I[A'(K,s)]\F s nA K , s } < 2{(6K) 1 / 2 + K}exp{±\ (l + e K )t}. 
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Proof: Note that, from Theorem 12.11 (1). given any e > 0, 

H' £ := supe- Ao(1+e) *||iJ(t)|| < ooa.s. . (2.22) 

This in turn implies that 

\h(t) - 1| = \Z{H r{d) {t)) - H r{d} {t)\ < sup \Z{u) - u\ for all t > 0, 

0<u<H' £ exp{A (l+e)t} 

and Part 1 follows from the law of the iterated logarithm for the Poisson process. In 
similar fashion, from (I2.16p . we have 

E{\h(t) - h(s)\I[A'(K, s)] | T s n A KiS } 

< E< sup \Z{u)—u\ JF s C\Ak,s 

[o<w<||i?(s)|| exp{A (l+£K)(t-s)} 

< 2(^) 1 / 2 exp{|A (l + ^)t}, (2.23) 
for Z{u) := Z{u + H r ^(s)) — Z(H r ^)(s)). Since also, on Ak, s , 

\h(s)-l\ < H rid) {s)^ 1+£ ^K^ 1 -^ < Kexp{lX s(l + e K )}, 
the proof is completed. ■ 
Recalling (12TT0|) . and writing W*(t) := e~ x ° t l T H(t) 1 we have 

dW* 



dt 



X e- Xot h{t), (2.24) 



and, in view of Corollary l2.3l (l). it follows that W*(oo) := lim^oo W*(t) exists and is finite 
a.s. The process W*(-), although directly motivated from the differential equations ( 12.101) . 
is not the usual choice for defining such a limit: the branching process martingale is 

r(d)-l 

W(t) := e~ Xot J2 H M = W*(t) + e- Xot {H (t) - H r(d) (t)} = W m (t) + e- Xat h(t), 

j=0 

this last by (12. lip . Thus lim 4 _^ 00 W / (t) = W*(oo) also, because of Corollary 12.31 (1). Note 
that, in similar fashion, (12. 24ft can also be written as 

^ = \ e- Xot {H (t)-H r(d) (t)}, (2.25) 
from which, by partial integration, it follows that 

POO 

W*(oo) = W*(0)+ / e- Xot {dH (t)-X H r{d) ^(t)dt}, 
Jo 

identifying W*(oo) as r(d)W, where W is the limiting random variable defined in |4], 
Theorem 4.1. 
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Similarly, integrating ( 12.1 Op . and noting that H(0) = because M/(0) = for all I > 1, 
we have ^ 

e -*> c rw*H(t) = [ X C r{d) e~ x ° c ^ u e r{d) h(u)du, 
Jo 

from which it follows that 

e- Xot H(t) = e- Xot [ X C r{d) e XoC ^ {t - u) e r{d) h(u)du. (2.26) 
Jo 

Since E/t(it) = 1 for all u > 0, it follows directly that 

e- Ao *||E#(t)|| < e- Ao * f X e Mt - u) \\e r{d) \\du < 1, (2.27) 

Jo 

so that e~ Xot H(t) is uniformly bounded in expectation. With some further effort, the long 
term behaviour of the vector e~ Xot H(t) can be related in detailed fashion to that of W*(t). 

Theorem 2.2 As t — > oo, W(t) := r(d)e~ Xat H{t) — > W*(oo)l a.s. Furthermore, for any 
< s < t, 

(1) \\W(t)-W(s)\\ = 0(exp{-i(l-£)A oS }) a.s. 
for any e > 0, if r(d) < 6 ; and 

(2) \\W(t) - W(s)\\ = 0(exp{-(l - cos(27r/r(rf)))A s}) a.s. 
if r{d) > 7. Finally, for A'(K, s) as defined in $2.161) , 

(3) E{||W(t)-W(s)||/[A'(A:,s)]| J s n4,J = O(#exp{-/3 r(d) A s}), 

w/iere & = \{l - e K ) if d < 6, and & = (1 - cos(2vr/rf)) if d > 7 and K 1 ^ > 
5/{2cos(2tt/7) - 1}. 

Proof: Since the eigenvalues of C r are the r-th roots Uj := e 2n U- 1 ) l / r of unity, 1 < j < r, 
and the corresponding eigenvectors e^ := r -1 / 2 ^^ 1 , u~ 2 } . . . , uj^ r ~ x \ 1) are (complex) 

orthonormal, we can write e~ Xot H(t) = Y7j=l ^j(t)e^ with faj(t) := e~ Xat {e^} T H(t), and 
use (12.261) to compute the hj(t). 

Taking j = 1 first, pre- multiplying by r(d) 1 ^ 2 e t ^ = 1 gives 

r(d) 1/2 /n(t) = e- Ao * f X e Xo{t - u h T e r(d) h(u)du = f X e- X ° u h(u) du . 

Jo Jo 

This last converges a.s. to W*(oo) in view of (12.241) . and hence h\{i) — > r(d)~ 1 ^ 2 W*(oo) 
a.s. as t — > oo. Indeed, we have a little more: for any s <t, from Corollary 12.31 (1). 

r(d) 1/2 |/ii(*)-/ii(s)| < / A e- Aou |/i(M)|rfM < — — exp{-±(l -e)X s} a.s. , (2.28) 
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for any e > 0. 

For 2 < j < r(d), we have 

r(d) 1/2 /i,(t) = e- Xot [ Xoojje^^h^du, 

Jo 

giving 

r(c/) 1/2 |^(t)| < e" Aot f X e Xop ^ u) \h{u)\du, 

Jo 

where pj := cos{27r(j — l)/r^)}- The estimates given in Parts 1 and 2 now follow directly 
from Corollary 12.31 (1). and 

\ime- Xot H(t) = r(d) -1 W,,(oo)l a.s. (2.29) 

t— >oo 

is immediate. The bound given in Part 3 follows in a similar way, but using Corol- 
lary 0(2) in place of Corollary 0(1). ■ 



Note, in particular, that the distribution of W*(oo) is the same, for given value of r(d), 
irrespective of the value of Aq, since W(u/Xo) = r(d)e~ u H (u/Xq) = r(d)e~ u H(u), with H 
as defined in Remark 12.11 

For use in Section El we define 

:= E{e- ew * ((Xl) | M o (0) = 1, M,(0) = 0, j > 1}. (2.30) 



The function 0, being the Laplace transform of a branching process limit random variable, 
can as usual be expressed as the solution of an implicit equation. This is based on the 
observation that, because of the branching property, 

W*(oo) = d ^e-^'W®, 

where (W* , j > 1) are independent copies of W*(oo), and (tj, j > 1) are the event times 
in a Poisson process on R + . For the gossip process in dimension d, the Poisson process 
has intensity 

pv{K)u d = X d+l u d /d\, u>0. 
This implies that the Laplace transform (f) satisfies 

exp(-/ 1{1 -0(9e- XoU )}X d o +1 u d du 



o 



OO 1 

= ex P|-y -{1 - (f)(9e- x )}x d dx j , (2.31) 

with -0'(O) = E{lV,(oo)} |M (0) = 1,1^(0) = 0, j > 1} = 1. If we define h(t) := 
1 — 0(e 4 ), then /i satisfies the equation 

h{t) = 1-expj-jf ^-/i(t - x) dxj , (2.32) 
12 



with \im t ^ ) ._ 00 e~ t h'(t) = 1; in the case d — 2, (I2.32p is the equation that appears in 
Chatterjee & Durrett ([5], Lemma 1.1). Note that the functions <p = <pd and h = hd thus 
depend only on the dimension d, and not on the choice of neighbourhoods, and hence 
that this is true also of the distribution of W*(oo). 

It follows from (I2.3ip that the lower tail of W*(oo) can be easily bounded: 

F{W*(oo)<w] < e(f)(l/w) = eexpj-^ -^{1 - <p(e~ x /w)}x d dx 

< eexpl- —{l-<p(l)}x d dx\ 

< e(exp{-c(log(l/ W )) d+1 }), (2.33) 

for c := {1 — 0(1)}/ (d+ 1)!. Thus P[W*(oo) < w] goes to zero faster than any power of w 
as w — > for all d > 1. For the upper tail, 

P[W*(oo) >w}< l/w (2.34) 

is immediate from Markov's inequality. 

Our final result of the section is a lower bound on the growth of H, needed later to 
show that F[A£' a ] is large enough, if s is not too small. 

Lemma 2.4 Let denote the time taken for H r ^d) to first reach a given level K . Then 
there is a constant Cc ■= (r(rf)!) 1 / r W/ 1 og(6/5) such that 

for all R>0. 

Proof: Define t := Xq 1 (r(d)\) 1 ^ d \ so that, starting X* with a single particle at 0, we 
have H r (d)(to) > 1. Then it is immediate that the process M is stochastically bounded 
below by a process M , where (M (nt ), n > 1) is a Galton-Watson process with offspring 
distribution pi — (1 — p 2 ) = e -1 , and, for nt Q < t < {n + l)i , M (t) = M (nt ). 
Furthermore, for all n > 1, H r (d){nto) > M (nto). Now, since Bi (m, i){[m/2,m]} > 1/2 
for all m > 1, a calculation shows that 

E{l/M ((n + l)t )\M (nt )=m} < m^f + f), 

from which it follows that ((6/5) n /Mo(nto), n > 0) is a non-negative supermartingale, 
with initial value 1. So, defining uk '■= min{n: M (nt ) > K}, and noting also that 
^o(^k^o) — 3K/2, it follows from the optional stopping theorem that 

E{(6/5)" K } < 3K/2. 

It thus follows that F[v K > i?/log(6/5)] < (3/2)^"^, for any R > 0, and hence that 

P[T^ (d) > t R/ log(6/5)] < (3/2)i^e- B 

also. ■ 
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2.2 Intersection asymptotics 

The branching process X* gives a useful approximation to Y as long as it has only few 
pairs of self intersecting islands. Thus our asymptotics for the branching process are of 
most interest over the time scale until self-intersections in X* become plentiful. To see 
when this is, we use (11.21) and (jl.3p . describing the intersection probabilities of randomly 
placed islands, to derive formulae for the mean number of pairs of self intersecting islands 
of X* at a given time t, and the number of islands of one process X{ that intersect islands 
of another, independent process X^. 

We begin by supposing that c g = 0. In this case, the number N(t) of self intersecting 
pairs of islands in a branching process at time t, conditional on having M (t) = n contacts 
taking place at times r*, . . . , r* < t, and with Tq = as usual, has mean given by 

E{N(t)\M (t) = n,r^...,T*} = n(n + V,t,t-Tl,...,t-r*) 

n— 1 n d / r\ 

i=0 j=i+l 1=0 ^ ' 

= §L-M/C) (f) WW^W - M d (t)}- (2-35) 
i=o ^ ' 

Similarly, again if c g = 0, for two independent branching processes, one that has developed 
to time t and has Mo{t) = m contacts taking place at times r*, . . . , < t, and the 
other that has developed to time s and has M (t) = n contacts taking place at times 
a*, . . . , cr* < s, the conditional mean of the number N(t, s) of intersecting pairs is given 
by 

E{N(t, s) | M (t) = m, t*, Mo(s) = n, a*,, . . . , <} 



n'(m + 1, n + 1; t, t — , . . . , t — r^; s, s 



a 



i ■ 



i=0 j=0 1=0 ^ ' 

= L^vW^f^MMM^a), (2.36) 
i=o ^ ' 

where Mi(s) := ^™ =0 (s — cr*) 1 . Thus the quantities Mj(-) of the previous section are 
exactly the quantities needed for making such computations. If c g > 0, we have 

\E{N(t) | M (t) = n,rt, -..,<}- fi(n + 1; t, t - r{, . . . , * - r*)| 



< ^Jc g {(2t)^(/C)/L} 1+ ^/ d ; 

E{iV(t, s) | M (t) = m, r*, . . . , r* , M ( s ) = n, a*, ...,<} (2.37) 

- //(m + 1, n + 1; t, £ - r*, . . . , t - r* m ; s, s - o\, . . . , s - a* 

< mnc g {{2t) d v{JC)/L} 1+ ^ /d . 
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Now, from Theorem 12.2} we have 

M^M^t) ~ l\(d-l)\\ d r(dy 2 e 2Xot , 
so that the mean number of self-intersections is small up until times t at which 

e 2X Q t x A ;= L x d Q / v ()C), (2.38) 

or, equivalently, for 

t = t A)X := -L{l gA + x}. (2.39) 

As can be seen from (I2.10p . Ag 1 is the time scale for the growth of H(-), or equivalently, 
in view of ( 12. 7\i and ( 12. 8p . the time scale on which the first new contact occurs. Corre- 
spondingly, v(IC)\Q d represents the scale for the size of the initial island at the time when 
the first new contact occurs. For the sort of asymptotics to be discussed here, with many 
islands involved in the covering of C, it is thus natural to require that v(IC)X d < L; 
hence we shall think of the dimensionless quantity A as being large. 



3 The deterministic phase 
3.1 Outline 

In this section, we show that the development of the proportion of the volume 

of C that is covered at time t grows more or less deterministically, once the initial stages 
have passed. We no longer make the simplifying assumption that c g = 0. As noted above, 
up to times of the form |AQ 1 {logA + x}, for x fixed, there are few self-intersections in 
the branching process X*, so that calculations made using the branching process can be 
expected to give close to the right answers for the small world and gossip processes as 
well. 

This fact was exploited for small world models in [4], when approximating the dis- 
tribution of the distance dsw(P, P') between two randomly chosen points P and P' in 
the small world graph. The key observation is that dsw{P, P') > 2t exactly when the 
sets Yp(t) and Yp>(t) are disjoint, and the discussion above indicates that the sets Yp(t) 
and Yp/(t) can be replaced with little error, for the calculation of probabilities, by the 
sets Yp(t) and Yp,(t) generated by independent branching processes X* P and X* pl — as 
it happens, provided only that intersections between islands J and J' of the branching 
processes are not counted when either J C J' or J' C J, since such a constellation cannot 
occur in the small world process. Thus, for P an independent uniform point of C, we 
have 

Fp [d sw (Po,P)<2t] = p[Y Po (t)nYp(t)^0] ~ P[y; o (t)n*y;(t)^0], (3.1) 

for independent Xp Q and Xp, where the notation n* is used to denote this special mode 
of intersection. The complement of the latter probability was approximated, for small 
world processes, in jl], Theorem 4.2. The connection with the current problem is that 

¥ Po [d G (P ,P)<2t} = L- 1 EV Po (2t), (3.2) 



15 



and indeed that 

L- l ¥.{V Po (2t)\F s } = F Po [d G (Po,P)<2t\F s } ~ F[Y* o (t)n*Y*(t)^(H\T s } =: ir(s,t), 

(3.3) 

where the asymptotic equivalence can be expected much as for (13. ip . Here d G (P,Q) < 
t -<=>- Q G Y P (t) describes the shortest distance between two points P and Q in the 
gossip process. Our aim is to demonstrate that the quantity L~ 1 Vp (2t) stays close to 
its expectation L -1 E{Vp (2t) \J- S }, once a certain time s has elapsed, showing that the 
volume process develops in almost deterministic fashion from s onwards. For this, it is 
enough to show that the conditional variance Var {L~ 1 Vp Q (2t) | Tg) becomes small with s. 
Now the expectation L~ 1 K{Vp (2t) \ J- s } can be approximated, as above, by n(s,t). For 
the mean square, we simply note that, for t > s, 

L- 2 E{[V Po (2t)] 2 \ F s } = F Po [{d G (P ,P)<2t}n{d G (P ,P')<2t}\T s ] 

~ p[{F; o (t)n*y;(t)^0}n{r; o (t)n*r;,(t)^0}|j-f o ], (3.4) 

using three independent processes Xp o , X* p and Xp, , where P and P' are two independent 
uniform points of C, and J 7 ^ denotes the history of X Pq . Since, for s large enough, the 
development of the branching process X Pq after s is almost deterministic, the statistics of 
the set Yp (t) are already almost determined at time s. Hence the probability in (13. 4p is 
close to 

{P[y; o (t)n*r;(t)^0|j-f°]} 2 = : 7r 2 ( s ,t), 

asymptotically equivalent to the square of the conditional mean. In consequence, the 
conditional variance is small, as required. The remainder of the section consists of making 
this approximate argument precise, in the context of gossip processes; note that, for these 
processes, the special intersection D* is not required. For small world processes, the 
corresponding estimates can be deduced using the methods of [1], Section 4. 

3.2 Constructions 

For the gossip process, the following construction of Y is useful. First, for any (P, t) £ 
C x [0, oo) let S;c(P, t) C C x [0, oo) denote the set whose C-section ctc,u{Sk(P, t)} at u is 
/C(P, u — t) if u > t and otherwise. Let n denote a marked Poisson process on C x [0, oo) 
with constant intensity p, and with marks independently and uniformly distributed in C. 
Take P G C to be the initial point; set So := S/c(Po, 0). Then define Si to be the set 

5*1 := S [J S/c(Py, Tij), 

where the points of n in So occur at locations and times (Qij,Tij) and have marks Py, 
j > 1. Then recursively, for I > 2, define 

Si := Si-i S)c(Pij,Tij), 

where the points of n in Si-i \ S1-2 occur at locations and times (Qij,Tij) and have 
marks Py, j > 1. Then we can define 

Y(t) := acA S ao}, 0<t<T, (3.5) 
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where := U,> S/. Clearly, D {C x [0, T]} = SiD{C x [0, T]} for some I > 0, since n 
has a.s. only finitely many points in C x [0, T\. 

Thinking of the points (Qij,Tij) with marks Pij as the Z-th generation descendants of 
the initial individual Pq, we note that the process Y can be identified with the Markov 
branching process X* until the first time t at which there is a non-empty intersection 
between two of the sets Sjc(Pij,i), for I > and j > 1 such that Ty < t; this is because 
of the independence of Poisson realizations on disjoint sets. Thus the construction of Y 
also yields a part of a coupled X*. 

As an alternative, the construction of II and Y can be replicated starting from the 
branching process X*. Start with P , and, writing So := Sjc(Po,0), assign Poisson points 
(Qij,^ij), j > 1, to So fl {C x [0, T]} with intensity p, and with marks independently 
and uniformly distributed on C; denote these by Jy := (Qij, fy; Pij), j > 1- These are 
the first generation descendants of the individual at Pq that are born before T. Repeat 
the process recursively, at the Z-th. step, / > 2, assigning Poisson points and marks to 
each fl {C x [0,T]}, where Si-ij := Sfc(Pi-i ) j,?i-.i ) j), yielding the Z-th generation 

descendants ly := (Qij, %; -P/j), j > 1, of the individual at P that are born before T; as 
before, write Si := 5j_x U {Uj->nSy}. There are in general more of these descendants than 
there are Poisson points and marks in fl {C x [0,T]}, and the labelling is typically 
different. To recover the Poisson process II and its marks, first identify the set {hj, j > 1} 
with {(Qij, Tij] Pij), j > 1}, the first generation descendants of the individual at P in the 
gossip process. Then label each with Gij G {0, 1}, according as whether the point is to 
be treated as real or a ghost; Gij = 1 if (Pij, Ty) G 5*0, and, if there is any j with Gij = 1, 
the union in (13. 5p is not disjoint beyond Tij. 

For the remaining construction, the descendants J/ J;r of individual I\j in [0, T] are the 
Poisson points and marks from Sij fl {C x [0, T}} chosen above. If Gij = 1, label them all 
with Gij ]r = 1 — as descendants of ghosts, they are themselves ghosts. If Gy = 0, label 
those of the ly ;r . for which 



with Gij ;r = 1 also, since they lie in a part of C x [0, T] that has already been covered, 
and therefore do not belong to II. Do the same if 



since such points do not generate new descendants in Y, because then Sic(Pij-r,Tij-r) 
is contained in the already covered region of C x [0, oo), even though here the points 
(Qij-r, Tij- r ) do represent points of II. The remaining points determine the points of IT in 






carry the label Gij ;r = and are in C x [0, T]. The set of 



points 



(Qlj;r,Tlj- r ), j > 1,T > 1 : (Qlj- r , fy ;r ) ^ [J Slj> 

j'=l 
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together with their marks fy ;r , can now be identified with the set of points and marks 
{(Qi+i,ji T i+i,j] Pi+i,j)i j > I}; after a suitable re-indexation (and with the labels Gi + ij 
correspondingly defined), recovering the gossip process Y. 

Note that the same construction can also be used starting from two independent 
branching processes X^* and X^* with initial points P an d Pq\ provided that, in 
each generation, the descendants of both individuals are listed together. The resulting 
gossip process describes the informed regions at each time t, when the information spreads 
from two initial sources at Pq and P^ ' . Note that the order in which individuals appear 
in the list can influence the realization of Y that is obtained, so that, for instance, the 
set Yp(i) (t) obtained from X*- 1 -** alone may not be identical to the subset of points within 

distance t of P m the set Y (i) D (2)(t) obtained from X^* and X^* together. However, 

' 

it is shown in the proofs below that the differences, which only occur as a result of 
sets in the branching processes overlapping, are not significant for the ranges of t under 
consideration here. 

For the calculations to come, we also need an analogue of (13. 2p for gossip processes. 
Here, the argument is a little simpler than for small world processes. A point P has been 
informed from Pq at time 2t exactly when the set Yp (t) C C of points informed from Pq by 
time t intersects the set of points Yp(t; 2t) C C from which the information will reach P 
by time 2t if it has reached Yp(t; 2t) by time t. Now the set Yp(t; 2t) is determined by 
the points (rj, Pj) of n in C x (t, 2t], together with their associated marks Q*, and can be 
constructed from them in exactly the same way as Yp (t) is constructed from the points 
(Tj,Pj) of n in C x (0, t], together with their associated marks Q®, except that time is 
run backwards from 2t to t, and the roles of Pj and are swapped. Thus, and because 
the neighbourhoods /C were so chosen that P G JC(Q, t) exactly when Q 6 JC(P, t), the set 
Yp(t;2t) has the same distribution as Yp(t), and Yp (t) and Yp(t;2t) are independent. 
Thus we have an analogue of (I3.2p for gossip processes: for any s < t, 

L- 1 E{V Po (2t)| jr s } = W>p [Yp (t)nYp(t;2t)^®\F s ]; (3.6) 
L- 2 E{[V Po (2t)} 2 | F s } = ¥ Po [{Y Po (t) n Y P (t; 2t) ± 0} n {Y Po (t) n Y P ,(t; 2t) ^ 0} | T B ], 

where P and P' are independent uniform points of C, and the argument is now primarily 
aimed at justifying the replacement of Y Po (t), Yp(t;2t) and Yp/(t;2t) in the formulae 
by independent copies of Yq{t) with appropriate choices of Q, so that computations can 
conveniently be made. 

3.3 Calculations 

As remarked at the start of Section [21 it is useful to be able to bound X* above and 
below by branching processes having c g = 0, to which we can then apply the results of 
that section. We do this by constructing processes X~ and X + , for times t > s, with the 
same initial conditions as X* at s, and using the same underlying Poisson process Z; the 
time s and the initial conditions can be freely chosen. For X*, the quantity v{K)M^{u) 
has to be replaced by V*(u) in the argument of Z in (12.41) . and, by (II. ip . we have 

v{K)Md{u){\ — ?7a} < V*(u) < v()C)M d (u){l + VA}, 
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where 

/31ogAV* 

VA := ° g \2A^ J ' ( } 

uniformly in < u < (3/2)Aq 1 log A. Hence we can define X + as in Section [2] by using 
p + := p{l + T]\} as the contact rate per unit volume, and X~ with p~ := p{l — ?7 A }, in 
which case 

Mf(t) < M t (t) < M+(t) for all < t < (3/2)A 1 log A and < I < d+ 1. (3.8) 

We shall use X~ and X + extensively to bound quantities associated with X*, and write 

A± := X {l±r ]A } 1 / d (3.9) 

for the corresponding growth exponents; for convenience, we shall assume henceforth 
that A is large enough that 9rj^\ogA < 1. 

We now continue with the following Poisson approximation result (see jl], Theo- 
rem 3.1), which can be simply proved using the Stein-Chen method. 

Lemma 3.1 Let n K-islands of radii ti,...,t n < (3/2)AQ 1 logA have centres indepen- 
dently and uniformly distributed on C , and let N n denote the number of pairs of them 
that intersect. Then 

d TV (£(N n ),Po(EN n )) < 4np+ 

where, recalling (GUP, 

Pt M{ 1 + C9 (^'} = M {Wa} ,3,0, 

is an upper bound for the probability that two independently positioned JC-islands of ra- 
dius at most (3/2)AQ 1 logA intersect. Similarly, for two independent collections of K- 
islands, one with m and one with n islands, having radii £i, . . .,i m < (3/2)Ao 1 logA and 
Ui, . . . , u n < (3/2) Aq 1 log A respectively, the number N mn of intersecting pairs satisfies 

d TV (£(N mn ),Po(EN mn )) < 2{m + n)p+. 

Thus, when finding the probability of there being an intersection between the islands 
of X Pi (t) and Xp 2 (t), where X* Pi and Xp 2 are independent, this Poisson approximation 
offers an approach. It is exploited in the following result, in which AfW and denote 
the quantities f)2.2p derived from X Pi (t) and Xp 2 (t). 

Lemma 3.2 Define 

P*[m (1 W 2) ] := P[^(t)nr; 2 (t) ^0|M«(t) = m«,M^(t) = m®\. 

Then, for t < (3/2)Aq 1 log A, we have 

P *[ m M m (2)] _ (i _ ex P {-]T (^mfm^L-^lC)}) 
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In particular, fort = t^ x := |A 1 {log A + x}, as in \2. 39^1 , forx < |logA and for K > 1, 
this gives 



P*[MM(t),MW(t)} - (l - exp{- £ ( d )M^\t)M^(t)X d A^} S 

< C |irA-5+! £ A'{ log A} d e x/2 + K 2 A-^ /d+3£K/2 {\og A} d +^j , (3.11) 
for suitable choice of C, except on an event of probability at most 2c a e~ Rl 3//5 . 



Proof: The first part of the lemma is an immediate consequence of Lemma I3.1[ together 
with (I2.37p . For the second part, we use (13~8|) to bound (t) and M (2) '(*) above by 
M (1)+ (t) and M^ )+ (t), and then Theorem [0(1) to show that, for I = 1,2, 

M (0 (t) < C a Ke A o' (1+£K) < eC a Ke Xot{1+£K) 

for all < t < (3/2)A 1 logA, except on a set of probability at most c a e ' , since, 
for t in this range, 

A+t < Aot + (3/2)^ A logA(l + £ X ) < A t + 1, (3.12) 

if 9t]a log A < 1 (since also Ek < 5 for K > 1). ■ 

We now aim to show that, for t = t^ )X and s < t, the conditional probabilities of actual 
interest for (13. 6p . 

Pl {s,t) := F Po [Y Po {t)nY P {t;2t) J"J (3.13) 

and 

p 2 (s,t) := Pp [{F Po (t)nFp(t;2t) ^0}n{F Po (t)nF P ,(t;2t) ^0}| J",], (3.14) 

are close to probabilities 

Pl(s,t) := E Po {P*[M^(t),M^(t)]\T s }; 

p*(s,t) := EpJP^MW^^^^P^MW^),]^^)^)]! F 8 }, 

that we can approximate using Lemma 13.21 Here, is used to denote the quan- 

tities (12. 2D for X Po , coupled as above with Y Po , and J 7 , to denote its history up to s. 

and are related to Y P (t;2t) and Y P /(t; 2t) in similar fashion, through branch- 
ing processes Xp and Xp,, which are independent of each other and of X Pq . Now the 
union Yp o (t) of the islands of Xp Q (t) contains Y Po (t), and the corresponding statement 
is true for X P (t) and Y P (t]2t) and for Xp,{t) and Y P /(t;2t). From this, it follows that 
P*( s ,t) > pi(s,t) and that p^s.t) > p 2 (s,t). We thus need only to show that the differ- 
ences pl(s, t) —pi(s, t) and P2(s,t) — p 2 (s,t) cannot be too large. This we establish on the 
event Ak s £ J~ a , defined as in (l2~PTj) with Mf \t)Xi/i\ for H^t), for a suitable choice of 
K = K(A). 
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The differences between Pi(s, t) and pi(s,t), I = 1,2, arise from events on which an 
island of a branching process X* makes an intersection that is not an intersection in 
the corresponding ^-process. Such an event can only occur if the island of X* is a 
ghost, or if the intersection occurs at a part of an island that overlaps another island, 
so that an intersection may have been counted twice using X*. Thus it will be enough 
to show that the probability of there being a ghost or an overlapped island in X Po (t) 
that intersects Y P {t;2t) U Y P /(t;2t) is small, and that the same is true for ghosts and 
overlapped islands of X* ppl (t) := X* p {t) U X* pl (t) intersecting Y Po (t). 

Lemma 3.3 Define K(A) := (40 log A) 3 . Then, on A K ( A ^ S , for t A)X := ^{logA + x} as 
in ( TO^) . and with \x\ < (1/6) log A, 

< Pt{s,t)-pi(s,t) = 0{e x {l + e x )A- 5 / 12 {logA} M+12 + A- 1 / 6 {logA} d+6 ), 1 = 1,2. 

Proof: To make the necessary estimates, we begin by coupling X Po in t > s to upper and 
lower processes X + and X~ as in (13.81) . starting from the same state at time s. Define 
the event A' + (K(A),s) as in (12.161) . but for the process X + , and so with Aq for A and 
with H + for H. Then, on the event A' + (K(A),s) n A^ jS , M+(t AjX ) < e6K(A)A 1 / 2 e x / 2 , 
and, in view of Theorem 12. 1^ A', (K(A), s) has conditional probability at least 1 — c 2 A~ 8 
on Ajc(a),s- Hence, the mean number of intersecting pairs of islands of X P[) at £a,x satisfies 

E{N(t A , x )I[A' + (K(A),s)}\ T s nA K{A) , s } < l{CK(A)A 1 / 2 e x / 2 } 2 pt < C e x {\og A} d+6 , 

(3.16) 

for some constants C,C, and also, from ( 13.101) . 

4E{Mo(*a,x)/K(#(A), s)} | T s n Ak{A), s }p1 < C"e^ 2 A^ 2 {\ogA} d+3 . (3.17) 

Thus, taking x = —(1/6) log A, the conditional probability of any pair of islands of X Po 
intersecting before time (5/12)Aq 1 log A is at most of order 0(A _1 / 6 {log A} d+6 ) on Akia),si 
from (13.161) . Also, if x < |logA, the number of such intersecting pairs before time t\ tX 
exceeds log A{1 Ve x {log A} d+6 } with conditional probability of order 0(A _1 / 4 {log A} d+3 ), 
in view of (13.161) . (I3.17p . Lemma |3.1[ and the Chernoff inequalities for the Poisson distri- 
bution. 

Thus, on Ak(a),si except on an event of probability of order 0(A -1 / 6 {log A} d+6 ), all 
intersections of pairs of islands of X Pq occur after time (5/12)Aq 1 log A, and, once more 
in view of Theorem 12. II applied to the dominating branching process X + , they each give 
rise to at most C(log A)^ d+1 ^ v3 ^A 1 / 12 e a: / 2 ghosts, except on an event of probability at 
most c a A~ 7 , since each such intersection at worst gives rise to a ghost branching pro- 
cess starting with a single island of radius (3/2)Aq 1 log A. Thus, on A K (a) >s , except on 
an event of conditional probability of order 0(A _1 / 6 {log A} d+6 ), there are at most of or- 
der 0{A l / l2 {\ogA} 2d+9 e x / 2 {l + e x )) islands in X* Po (t) whose intersections with Y P (t; 2t) U 
Y P i{t;2t) should be discounted, and their radii cannot exceed (3/2) A^ 1 log A. Further- 
more, the number of islands in X P (t) and X P ,(t) together is at most 

C a K(A)e x o^ 1+£ ^ = 0(A 1/2 {logA}V /2 ), 
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except on an event of probability at most c a A~ 8 , by Theorem 12.1( 1). Thus, on Ak(A),si 
the mean number of intersections that should be neglected is, off the exceptional events, 
at most of order 

O (A^^logAj^+V/^l + e^A^llogAjV'/VA) = 0(e x (l + e x )A- 5/12 {\ogA} 3d+12 ). 

This expectation, together with the complementary expectation from overlapping islands 
and ghosts in X Ppi , bounds that part of the differences p*(s,t) — pi(s,t), 1 — 1,2, aris- 
ing off the exceptional events, and the exceptional events together have probability of 
order 0(A _1 / 6 {log A} d+6 ). The argument for this second expectation is the same, except 
that there is no conditioning, making it equivalent to the previous argument with s — 0, 
and with A(K(A), 0) being automatically satisfied because the inital state consists of just 
two singletons. ■ 



Lemma 3.4 For M' 1 ' , and as in Lemma \3~3\ and for x < |logA ; we have 

d 

E (s) ' 



> ( ex p{ - E ( d ) M ? ] (tA,M^ } 



j\{d-j)\wf\s)e 



and 



< Ce x {e~^ d+lXoS (logA) 3 + (log Af + ^A~^ /d }, 
E(-) { exp { - E M f ] (^) [M^j (t^) + (t A , x )\ X d A- 1 } 

^ i'=n V' 



(3.18) 



(3.19) 



< Ce a; (logA) 6 {e- /3d + lAoS + (log A) 1+7 *A^ /d } J 

/or any s < t\ tX , where E( s ) denotes expectation conditional on J 7 , fl Ak(A),.s; H^- (s) := 
(d+ l)e _A ° s M- 1 ' ) (s)A^/j!, zs as in A2.30\) . and C,C are suitable constants. 

Proof: We give the proof of H3 . 1 9D : that of (13.181) is simpler. For t > s, we bound 
above and below by and M®~, 1 < / < 3, as in the previous lemma. We then 

observe that, for I = 2,3, 

Mf\t K , x )M ( i ) _ j {t K , x )\ d A- 1 

< Mj 1)+ (^)Mf+(t AiX )(A+) rf A" 1 

= (d+ l)- 2 j\(d- j^^W^it^W^t^A- 1 

< (d+l)- 2 j\(d-j)\e( 1+ ^ l °z A +^W^ 

< (rf + l)- 2 j1(rf-jO!e^ {1)+ (t A , 3; )^i J + (tA,,){l + 2er ?A logA}, 
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where r/ A is as in (13. 7p . We now use Theorem 12.21 (3) to deduce that 

E W -W^ + (s)\I[A' + (K(A),s)]) = 0{K{A)e~^^ s ) 

^ s) (\W^(t^)-W^(oo)\I[A' + (K(A),s)}) = 0{K{A)e~^% 



d 'd 



j=0 



E( S ) exp )Mf\t Kx )[Mf! 3 (t Kx ) + M^MA-i 



(3.20) 



/ = 2, 3, and Theorem 0(2) to show that F[A' + (K(A), s) | T s n A K{A)tS \ > 1 - c 2 A~ 8 . It 
thus follows, also using ( 12.27(1 . that 



> E W exp{-(rf+ l)- 2 d\ J2^ 1)+ (s)[W^ + (oo) + > + (oo)K 

^ 3=0 

- O {e x K(A) 2 {e-^ d +^ s + r] A log A} + A" 8 ) 
= {^(d+l)- 2 d\J2K\s)e*)} (3.21) 



j=0 



- 0(e*(logA) 6 [e-^ lAoa + (logA) 1+7fl A- 7fl/d ]), 

since Wj (s) = Wj (s) is .^-measurable, and since W* (oo) and W* (oc) are inde- 
pendent given J^s, and each has the distribution of W*(oo) as in f )2.30p . The upper bound 
is proved in analogous fashion. I 



Theorem 3.1 Choose s A '■= (oz/2)\q log A for some < a < 1/2. Let (3d be as defined in 
Theorem \2.2\ let Ak, s be as for Lemma l373[ and define t AjX as in j[2.39\) . Then, for any gos- 
sip process satisfying M.l\) , M.2\) and il.3\) , and for any 71 < 70 := min{|a/3d+i, \l g /d}, 
there exists a constant fc 7l < 00 such that 

Vax{L- 1 V Pb (2t A , a: ) I J SA nl K( A), SA } < KA-~<\ 

uniformly in \x\ < c^logA, where c v := A^g/d}. Furthermore, F[Ak(A),s a ] > 

1 — caA~ 72 for some 72 > and ca < 00. 

Proof: Starting from f)3.6p . collecting the results of Lemmas 13. 2 \ 13.31 and 13.44 and taking 
\x\ < I log A, we find after comparing the various errors that 

E{L- 1 V^(2t A , ai )|j; A nI A - (A ), aA } = G(W^(s A ),x) 

+ 0(e*{(log A) 3 A-°^+i/2 + (l og A) 4+7 «A- 7 « /d } + (logA) d+6+7 «A~ 7d/d ), 
= G(W w (s A ),x) + 0(e x (\ogA) 4+ ^A- 3c * + (log A) d+6+7 «A- 7 « /d ), (3.22) 

where 

d 

G{w,x) = l-(t) d {e x {d+l)- 2 d\^w^. (3.23) 

j=0 
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Similarly, 

¥.{[L- l V Po (2t Kx )f | F SA n A K{A)ySA } 

= {G(W {1) (s A ),x)} 2 + 0(e x (logA) 4+79 A- 3c " + (logA) d+6+7 *A- 7 ^). (3.24) 

Combining fj3.22j) and (I3.24p . and bounding the error term, the first statement of the 
theorem follows. 

To bound P[Ar-(a),s a ] f rom below, we start with Theorem 12.11 (1). Putting s = and 
then t = s A , and setting K = K(A) = (40 log A) 3 , it follows that 

F[e- x ^\\H(s A )\\ <C a K(A)e l / m ] > 1 - c a A~ 8 , 

so that P[A£L A) ] > 1 - c a A" 8 if 9 := C a e 1/80 . Then, from Lemma E31 

P[Ag{ A)iS J = F[H r{d) (s A )>K(A)] = F[r r ^ A) < s A ], 

bounded below by 1 — (3/2)K(A) exp{— saA /c c } = 1 — CA~ 7 for some 7 > 0, because of 
the choice of s A . Finally, from Lemma l2"T2"l (4). since 

ir(A)2( 1 -^(A))(2^(A)_i) x (bgA) 1 / 2 

exceeds 42 log 2 for all A sufficiently large, 

1 - n^\ A) , £K(A)>SA ] < c,ex P {-K(A)^^/28} 

< C ie xp{-C 2 (logA) 3 / 2 }, 
completing the proof. ■ 



This theorem is the basis for the main result of the paper, showing that the distribution 
of the path L~ l Vp Q (t) is concentrated close to its conditional mean. To complete the proof, 
we need first to have an expression for the conditional mean. 

Lemma 3.5 Uniformly in \x\ < c^logA, and with s A := (cn/2)Ag 1 log A for some < 
a < 1/2, 

EjL-Vp^AoHlogA + x})! J SA n^ (A) , SA } = Zi d (a; + logC , d + logiyi 1 )(s A )) + 0(A- 71 ), 

where ha is as in \2. 3ty) , Cd := (d + l)~ l d\, and 71 is as for Theorem IJ.il Also, if P is 
independently and uniformly chosen on C , the time Tp := inf{t > 0: P G Yp () (t)} satisfies 

¥[t p > Ao^logA + x)] = E(exp{-e^ d ^«(oo)^i 2 )(oo)}) + 0(A~ 73 ), 
for some 73 > ; where W*(oo) and Wl T (00) are independent. 

Proof: By arguing as for Lemma 13.41 using Theorem 12.2( 3). it follows that, uniformly 
for x as in Theorem 13.11 

E{L- l V Po (2t AjX ) I F SA n A K{A)jSA } 
= E{G(W^(s A ), x) I F SA n A K{A)tSA } + 0(A- 71 ) 
= E{G(W^(s A )l, x) I F SA n A K{A)iSA } + 0(A~ 71 ) 
= l-UC d e x W^ l \s A )) + 0{A-^) } (3.25) 
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where is as defined above, and the first part follows from the definition of hd follow- 
ing (I2.30p . The final result follows from taking the unconditional expectation in (13.251) . 
and applying Theorem 12.21 (3): 

l-E{L-V Po (2t A ,,)} = EWC/^Mll + Ofr'j+OfA^) 

= E{0 d (C d e^i 1 )(oo))} + O(A^ 3 ), 

for some 73 > 0. I 

Note that the form of the neighbourhoods /C only comes into the formulae of Lemma 1331 
through their volume v(JC), which is implicitly present in the time scaling by Ao in the 
definition of t A>x . 

We are now in a position to prove the pathwise approximation to L~ 1 Vp (t). Before 
doing so, we note that hd is the distribution function of a sum of independent random 
variables Z\ and Z 2 , where —Z\ has a standard Gumbel distribution, and Z 2 is distributed 
as — log W*(oo). This is because hd(x) := 1 — 0d(e x ) can be rewritten in alternative form 

as 

1-E je- eZ j = P[-Zi>-Z] = F[Z 1 + Z 2 <x], 
where Z := x + log W 7 ^ 00) is independent of Z\. 

Theorem 3.2 For any gossip process satisfying ( fi. 1}) . there exists a random variable U 
such that 

P[sup iL-Vp^Ao^logA + x}) - h d (x + logC, + U)\ > 4A" ai ] < cA^ 2 , 

X 

for some positive ai,a 2 and c < 00. e u has the distribution ofW*(pd), whose Laplace 
transform satisfies 112. 31}) , t A)X is as defined in 112. 39}) . h d is as in l}2. 32}) . and A is as 
defined in 112.38}) . 

Proof: Take s A := (a/2)Ao 1 logA for some < a < 1/2. It then follows from Theo- 
rem EUl Lemma [3.51 and Chebyshev's inequality that, for any a < j±/2, 

P[|£~ 1 V Po (Ao 1 {logA + x}) - h d (x + logCd + log WP(8 A ))\ > 2A~ a \ 7 SA n A KWt . A ] 

< k yi A 2a -^ (3.26) 

uniformly for \x\ < c„logA. 

Now, representing h using Z\ and Z 2 as above, we have, for x > 0, 

h d (-x) < P[-Zi > x/2] + P[logW*(oo) > x/2] < 2e- x/2 ; 
1 - h d {x) < F[-Z 1 < -x/2] + P[log W*(oo) < -x/2] < exp{-e* /2 } + e{x), 

where e(x) goes to zero super-exponentially fast, by (12. 33 p . Thus, for x\ := |c„logA, 
1 — hd(x\) = 0(A _1 ), and, for x~^ := — |c„logA, h d {x/ v ) < 2A~ C "/ 4 . Now choose an 

increasing sequence of m A := \A b ] points x^ between x~^ and x\, with x^' = x~^ and 
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X 



(tia) 



x A , in such a way that h d (x 



A ) 



helix®) ^ A~ fe for eacn i- On 



(A),s A ) 



from (BHD , 

-logW«(sA) < A sa = |«logA; logWW(s A ) < log(^(A)), (3.27) 

so that, choosing a < |c„, | logW* (s A )| < |c„logA for all A large enough. Thus we can 
take x = x A — log Wi (sa) in (I3.26[) for each < j < m A , obtaining 



P 



max |L-%,(Ao ^logA + xf - \ogW^\s A ) - \ogC d } - h d {xi 

0<J<m A 



> 2A~ 



F Sh n Ax(a),s A 



'"^ (3.28) 

< 2£; 7l A 2a+6 - 71 , 



from which it follows directly, because both Vp (t) and ^d(x) are non-decreasing in their 
arguments, that 



P 



sup |L-%(Ao MlogA + x}) - h d (x + \ogC d + log WW(s A ))\ 

> 2A~ a + A~ b 



*f(A),s A 



(3.29) 
< fcA 2a+6 ~ 7 \ 



for a suitable constant fc. 

The above estimate almost completes the theorem; it remains to replace W* (s A ) by 
iyi 1)+ (oo), where is as in Lemma 13.31 and then to remove the conditioning. For 

the former, we have 

E{\wPis A )-WP + ioo)\I[A' + iKiA),s A )} \ T SA nA KiA))SA } = 0(K(A) exp{-/3, +1 A+s A }), 



from Theorem CT3). together with P[A , + ( J fsT(A), s A ) | .F SA n ^(A), SA ] > 1 ~ c 2 A~ 8 , from 
Theorem 12.1( 2). Thus 

P[|^«(s A )-^i 1)+ (oo)| >A- b | J SA n^ (A) , SA ] = 0(A- b ), 

for any b' < \f3 d+ ia — b, and thus for any b' < 3c v — b, and hence also 

F[| log W«(«a) -log )+ (oo)| > A- fe / 2 1 J- SA nl x(A) , SA n{^i 1 ) + (oo) > A- 6 / 2 }] = 0(A- fo '). 

Since h d has density bounded by 1/e, the maximum of the density of the Gumbel distri- 
bution, it follows that 

P[sup |L~%(V{logA + x}) - ^(x + logC d + log^i 1)+ (oo))| 

- x 

> (2A~ a + A~ b + e-'A-o/ 2 ) | T SA n i^( A ), SA n {Wi 1)+ (oo) > A" b / 2 } 
< A;A 2a+fe - 71 + A;'A- 6 '. (3.30) 

But now, from fl2T33|) . P[W^ (1)+ (oo) < A~ b/2 ] = o(A~ u ) for any u > 0, and F[A K{A)tSA \ > 
1 — c^A -72 from Theorem 13. 1[ so that 



P 



sup |L-V Po (A HlogA + x}) - h d (x + hgC d + logiyi 1)+ (oo))| > 2A~ a + A" fe + e^A^ 2 

■ X 

< kA 2a+b -^ + C4A- 72 + c'A- y . (3.31) 
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Now take 2a = b = (c v A {7i/4}) and b' = c v to complete the proof, with ai = 6/2 and 
a 2 = min{7i/2,7 2 ,c 1 ,}. ■ 



For small worlds processes, the counterpart of Theorem 13.11 can be proved in entirely 
similar fashion. Lemma [3.51 is also correct, if is replaced by 

C d := dT 2 d\{{d +1)-1} = (d-l)!; (3.32) 

the change reflects both the difference in r(d) between gossip and small world processes, 
and the subtracted component arising from intersections forbidden in the small world 
context: see jl], Section 4 for more details. Note that the expression in |4], Theorem 4.2, 
appears different from that obtained here. It is, however, the same, with the transfor- 
mations r — y d, 2x — >■ x, ai — >■ v(IC)(f), W — > d -1 W*(oo), and noting that Lp = A/d\. 
For the limiting random variable W*(oo), the Poisson process of descendants of the first 
individual now has intensity 

pv(K,)du d - 1 = A^7(d-1)!, u>0, 

giving 

<P(0) = exp j- jT jj^y {l - cj ) {9e- x )}x d - 1 dx} , (3.33) 

the same equation as (12.311) for gossip processes, except that now d is replaced by d — 1; 
in view of (12.91) . this is not surprising. Thus Theorem 13.21 also holds, with h d replaced by 
h d _i and with Cd replaced by Ca- 
in the problem studied by Chatterjee & Durrett [5j, C is an iV x iV torus, so that 
d = 2 and L = N 2 , and /C(P, s) = B(P, s/V2n) is a Euclidean ball, so that v{K) = 1/2 
and c g = 0. They take p = N~ a for any a < 3, so that A = (d\pv(K)) l ^ d+ ^ = N~ a / 3 
and A = LX^/v(]C) = 2N 2 ^- a ^; also C d = 2/3 in Lemma [X3 Then the pathwise 
approximation given in Theorem 13.21 matches N~ 2 V(N a ^ 3 [2(l — a/3) logiV + x + log2]) 
with h{x + log(2/3) + U). This is seen to be the same as in Chatterjee & Durrett [5], 
noting that U has the same distribution as their logM, and that our choice of h(-) as a 
solution of (12.321) corresponds to their /z(- + log3), since they implicitly choose the solution 
of (12321) that has \im t ^ ^ e^h'(t) = 1/3. 

3.4 Complete coverage 

Theorem [372] and (13.251) . together with the fact that 2t\ x = Aq 1 {log A+x}, show that only 
a negligible fraction of C is covered at times much before Aq 1 log A, with C then becoming 
essentially covered in deterministic fashion on a time scale of order 0(1) around Aq 1 log A. 
Randomness is only visible in the time shift of — Aq 1 U in the origin of the transition from 
uncovered to covered. 

One can also, as in Chatterjee & Durrett [5], consider how long it takes until C is 
entirely covered. We show that complete coverage is achieved relatively soon after time 
Aq 1 log A, under the assumption that, for each s, C can be covered by n(s) islands of the 
form JC(P, s), where n(s) satisfies 

n(s) < c L/{v(IC)s d }, 0<s< L 1/d , (3.34) 
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for some c . 



Theorem 3.3 If Jff. 34\ ) is satisfied, then, except on a set of probability of order O (A 



-s\ 



for some 5 > 0, the whole of C is covered before time 

2 /721ogA\ 1/d 
T(A ' S)+ Aa^O ' 

where r(A, s) = A ^log A + 0((log A) 1/(d+1) )}. 

Proof: In a small world process, with probability 1— 0(A -5 ) for some 5 > 0, at least \Lp = 
\K/d\ shortcuts have been encountered by time r(A, s) := Aq 1 {log A + m — U — logC^}, 
where m denotes the median h~ 1 (l/2) and e u has the distribution of W*(oo), since then 
about half of C has been covered. In a gossip process, with similar probability, at least 
\A/d\ contacts have been made after time r(A, s) : = Aq 1 {log A + m — U — logC^} + Aq , 
since a volume of about L/2 has been generating contacts for a time interval of at least Aq . 
Write ip := d!/4, and note also that, from ( I2.33p . —U < (/clog A) 1 /^ 1 ) except on an event 
of probability of order 0(A -8 ), if k is chosen large enough. 

Cover C with islands of radius s = kX^ 1 ; with c as in (I3.34p this can be done using at 
most c L I '{v{K){k j X ) d } islands. Then, recalling f ll.ip . the probability that any of these 
islands contains none of Pi, ... , Py^a] is & t most 

c v(lC)A ( v()C)k d \ W coA ^ d/2 _ 



if A is large enough that c g (kA l l d )~ 1 9 < 1/2. On the complementary event, all of the 
islands, and thus all of C, are covered after an additional time of at most 2fcA ~ 1 . So take 
k = {(18/V>) log A} 1 ^, to make 

< -^tA- 8 ■ 

- 18 log A 



3.5 Manifolds with boundary 

The assumption that the manifold C is homogeneous simplifies the argument substantially. 
However, the adjustments needed if C is taken to be a 'reasonable' finite subset of a 
homogeneous manifold, such as a rectangle in R 2 or a spherical cap, are not great. The 
principal requirement is that most islands do not intersect the boundary dC. Let Cs '■ = 
{P G C: fC(P,8) fl dC 7^ 0} denote the ^-neighbourhood of the boundary of C, and 
assume that its volume is not too large when compared with that of C: 

L~ 1 \C 5 \ < c b 5(v{IC)/L) 1/d , 0<5<21ogA, (3.35) 

for some constant c&. For instance, with C a square of side yL in R 2 and with K.(P, s) 
a disc of radius s, d = 2 and v(fC) = n, and L -1 )^! < satisfies (I3.35P with 

Cfe = 4/a/7t. As is clear from the preceding argument, only times less than 5\ := 2Aq 1 log A 
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play a significant part, and the probability of an island with randomly chosen centre P G C 
intersecting dC before time 5\ is then at most 

c h S h {v{K)/Lfl i = 2c 6 A- 1/d logA (3.36) 

under the assumption (I3.35I) . 

In order to make the arguments of Section I3.3I work, it is enough to be able to bound 
the growth process above and below by branching processes with constant growth rates Xq 
and Aq , respectively, which are close enough to one another, and to have (II. 2p hold, with 
an error estimate similar to that in (II .3p . For the latter, if (jl.ip is satisfied, then (I3.36P 
implies that 

q L (t,u) = L^v(]C)(t + u) d (l + R L (t,u)), (3.37) 

where 

\R L (t,u)\ < c g {(t + u)(v(iq/L) l l d }^ +c b max{t,u}(v(IC)/L) 1 / d , (3.38) 

so that all that is needed is to replace the exponent 7 9 by 7 9 := min{7 9 , 1} and the 
constant c g by c g := c g + q, in (11.31) . when making intersection calculations. For the 
former, the upper bound Xq given in (13. 9 p still holds. A lower bound is obtained by 
neglecting any contacts to points of Cs A , and taking 

A " := A (1-L- 1 |C 5 |) > A {1 - 2c b A~ 1/d log A}, (3.39) 

by (I3.35P and (13.361) . Thus, once again, the previous arguments can be carried through, 
if, in the definition (13.71) . 7 9 is replaced by j g , and c g by c g + idcb/3. 

For the result corresponding to Theorem 13. 3[ it is also necessary to make the explicit 
assumption in connection with (I3.34p . which was previously guaranteed for all A such 
that c g {kh- 1/d )^ < 1/2, that each of the n(s) sets /C(P, s) used to cover C satisfies 

\K(P,s)nC\ > \s d v{K). (3.40) 

In the current context, this is no longer automatic, because part of a set JC(P, s) may 
lie outside C. The proof otherwise runs without any change. Note also that, under the 
extra assumptions (I3.35P and (13.401) of this section, the manifold C could also be allowed 
to consist of a number of disconnected components. In particular, the requirement (13.401) 
applied with s = kX^ 1 would prevent Theorem 13.31 from being justified if there were 
components that were too small, and this is to be expected, since it may take an extremely 
long time for a very small component of C to be hit by a sequence of randomly chosen 
points of C. 
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